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INTRODUCTION 


Free  space  laser  communications  provide  the  technology  which  will  revolutionize 
existing  information  networks  and  facilitate  the  creation  of  the  "information  super¬ 
highway".  The  information  networks  of  the  nearest  future  are  viewed  as  an  extensive 
system  of  satellites  connected  by  the  means  of  laser  communication  and  covering  the 
entire  planet.  Laser  communication  will  be  utilized  by  uploading  information  onto  the 
nearest  satellite,  creation  of  the  rational  satellite-to-satellite  communication  channel,  and 
downloading  information  to  the  required  location.  Utilization  of  laser  as  a 
communication  medium  offers  significant  benefits  over  traditional  RF  systems.  They 
include:  higher  bandwidth,  smaller  beam  divergence,  smaller  antenna  size,  weight,  and 
power.  Unfortunately,  the  smaller  beam  divergence  makes  the  system  susceptible  to 
performance  degradation  due  to  mechanical  vibrations.  The  feasibility  of  the  described 
information  network  depends  on  the  laser  positioning  technology  capable  of  offset 
adverse  effects  of  satellite  jitter.  This  work  addresses  a  small  but  important  part  of  laser 
communication  technology;  further  improvement  of  position  control  of  satellite-borne 
communication  lasers  by  the  application  of  advanced  control. 

Adaptive  control  techniques  can  significantly  improve  accuracy  and  dynamic 
performance  of  the  existing  pointing-acquisition-tracking  (PAT)  mechanisms  used  in  laser 
communication  systems.  Although  such  systems  utilize  the  most  advanced  materials  and 
technologies,  their  performance  is  adversely  affected  by  the  exposure  to  jitter  and  natural 
transients  (bending  modes)  of  moving  mechanical  parts.  The  main  motivation  for  this 
research  effort  is  in  the  fact  that  application  of  advanced  control  procedures, 
implemented  in  software,  can  further  improve  performance  characteristics  of  a  laser 
beam  steering  system  and  assure  its  high  performance  over  long  period  of  time.  The 
Model  Reference  (MR)  approach  has  been  developed  in  adaptive  control  as  a  technique, 
capable  of  "forcing"  a  control  system  to  follow  a  particular  performance  pattern 
represented  by  a  reference  model.  It  can  be  said  that  MR  approach  utilizes  the  feedback 
principle  for  "controlling  the  controller"  of  the  original  system.  As  such,  MR  results  in 
the  reduction  of  system  sensitivity  to  external  forcing  functions  (disturbances)  and 
parameter  perturbations  of  the  controlled  process.  While  MR  control  systems  exhibit 
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high  degree  of  robustness,  their  control  efforts  are  consistently  lower  than  that  of  robust 
control  systems  with  fixed  controllers.  At  the  same  time,  increased  amount  of 
calculations,  implementing  MR  control  laws,  does  not  constitute  a  difficulty  for  modem 
microprocessors. 

Originally,  the  MR  approach  was  suggested  for  the  modification  of  the  dynamic 
response  of  optical  mirrors.  It  allowed  to  extend  performance  limits  of  an  optical  mirror 
and  improve  its  dynamics  by  elimination  of  bending  modes.  MR  control  laws  can  be 
applied  to  reliable  but  inexpensive  microelectronic  devices,  often  just  by  the  modification 
of  software  of  already  existing  microprocessor  systems.  Therefore,  the  approach 
facilitates  application  of  relatively  simple,  reliable,  and  relatively  inexpensive  hardware 
with  performance  characteristics  significantly  extended  by  advanced  control  schemes. 

Effects  of  satellite  jitter  on  the  beam  positioning  system  are  quite  similar  to  the 
effects  of  bending  modes  of  the  optical  mirror.  The  MR  approach  implies  that  the 
response  of  the  existing  physical  system  is  "forced"  to  follow  the  response  of  a  simulation 
model  implemented  in  software.  The  model  is  subjected  to  the  same  commanded  input 
as  the  beam  positioning  system,  but  the  model  response  is  not  affected  by  jitter  (or  any 
other  undesirable  effects).  While  the  MR  scheme  eliminates  any  discrepancy  between 
the  responses  of  the  physical  system  and  the  model,  it  has  the  potential  for  jitter 
rejection. 

As  with  any  physical  system  subjected  to  long  term  effects  of  adverse  operational 
conditions  and  natural  wear,  performance  of  a  beam  positioning  system  is  expected  to 
deteriorate  with  time.  The  required  system  dynamics,  represented  by  the  simulation 
code  of  a  reference  model,  is  time-invariant.  Therefore,  the  required  performance 
characteristics  of  a  beam  positioning  system  can  be  assured  and  maintained  via 
application  of  the  MR  approach.  It  is  important  that  MR  technology  does  not  require 
any  additional  hardware  and  numerically  is  much  simpler  than  the  feedforward  adaptive 
jitter  compensation. 

An  investigation  of  MR  applications  for  reduction  of  jitter  effects  and  maintaining 
the  required  performance  characteristics  of  a  beam  positioning  system  over  long  periods 
of  system  operation  is  conducted.  Numerically  efficient  computational  schemes  of  MR 
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procedures  are  defined  and  tested,  both  by  simulation  and  experimentally. 

The  probability  of  burst  errors  in  an  intersatellite  lasercom  system  is  dependent, 
among  other  factors,  on  the  laser  power  and  laser  positioning  accuracy.  The  approach, 
presented  herein,  results  in  increased  accuracy  of  beam  positioning.  This  implies  that  the 
required  probability  of  communication  errors  can  be  achieved  with  lower  power  of  the 
laser  beam.  Operation  of  the  laser  source  at  a  "lower  than  nominal"  power  level  results 
in  the  improved  reliability  of  power  electronics.  In  addition,  although  solar  energy 
presents  an  inexhaustible  power  source,  reduced  power  consumption  can  be  beneficial 
for  extending  life  time  of  satellite  batteries.  Some  analyses  of  the  accuracy/power  trade¬ 
off,  leading  to  the  definition  of  a  control  law  governing  the  laser  power  as  a  function  of 
time-varying  beam  positioning  errors,  are  obtained.  Application  of  such  a  law  would 
result  in  the  optimal  use  of  the  available  power  sources,  and  increased  reliability  of 
power  electronics. 

A  further  investigation  and  application  of  the  suggested  control  law,  including 
analysis  of  the  distribution  laws  of  beam  intensity,  definition  of  probability  density 
function  required  for  the  analysis  of  burst  error  probability,  development  of  a  numerical 
procedure  and  computer  code,  is  recommended. 

The  experimental  part  of  the  proposed  studies  is  conducted  on  the  PAT  testbed 
of  US  Air  Force  Rome  Laboratory  with  the  participation  of  Dr.  Timothy  E.  Busch  [21]. 


1.  MATHEMATICAL  DESCRIPTION  AND  COMPUTER  SIMULATION 
OF  BEAM  STEERING  MIRRORS 


The  use  of  optical  frequencies  results  in  narrow  beams  and  less  beam  spreading  as 
the  beam  propagates  over  the  intervening  distances  between  satellites.  As  a  result,  a 
much  higher  degree  of  pointing  and  tracking  accuracy  are  required.  Beam  steering 
devices  are  used  to  provide  fine  steering  and  to  stabilize  the  line  of  sight  of  the 
transmitting  beam  from  platform  jitter.  The  devices  fall  into  the  two  categories  of  non¬ 
mechanical  and  mechanical  systems.  Mechanical  devices  are  mirrors  that  are  attached  to 
a  tip/tilt  platform  and  steer  the  beam  through  the  movment  of  the  mirror.  Two  types  of 
mirrors  are  commonly  in  use.  One  represents  the  state  of  the  art  in  high  bandwidth 
beam  steering  while  the  other  is  a  low  cost,  commercially  available  device  with 
substantially  poorer  performance.  Both  beam  steering  mechanical  systems  are  subjected 
to  mathematical  description  and  simulation  models. 

1.1  High  bandwidth  fine  steering  mirror 

This  beam  steering  system  includes  a  primary  surface  mirror,  two  axis  flexures,  a 
support  structure,  eight  actuators  and  four  position  sensors  [1,2].  It  is  intended  for 
stereeng  the  beam  along  the  azimuth  and  elevation  axes.  Moving  mechanical  parts  of 
the  mirror  are  balanced  so  that  the  mirror  assembly  would  not  impart  any  additional 
vibrations  to  the  space  vehicle. 

Mathematical  models  facilitate  the  investigation  of  various  control  system 
strategies  were  established  on  the  basis  of  a  series  of  laboratory  tests.  A  preliminary 
mathematical  model  of  the  mirror,  and  quadrant  detector  including  their  drive  circuitry, 
was  established  on  the  basis  of  input/  output  measurements.  This  model  was  then  used 
to  formulate  and  justify  by  computer  simulation  some  novel  control  techniques  [9,14,15]. 
However,  attempts  to  implement  the  controller  based  on  this  model  revealed  that  it  was 
unable  to  accurately  predict  the  systems  performance.  This  led  to  the  development  of  a 
much  more  detailed  model.  The  development  of  both  models  is  presented  below. 

A.  Simple  model  of  high  bandwidth  fine  steering  mirror.  During  the  laboratory 
tests,  voltage  signals  were  applied  to  the  azimuth  and  elevation  inputs  of  the  mirror  drive 
system.  The  responses  of  the  mirror  position  were  recorded  by  a  PC  with  the  sampling 
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period  of  25x10  ®  sec  connected  to  the  output  of  the  position  subsystem.  The  LAB  VIEW 
software  package  was  utilized  for  the  implementation  of  data  acquisition  system. 
Mathematical  model  of  the  mirror  dynamics  was  defined  as  follows: 


r  1 

r 

1  r 

|a(s) 

ID 

G*i!(s)  I  |v^{s) 

1 1 

|e{s) 

G„{S) 

1  1 

G„{s)  1  |v,(s) 

L  J 

L 

J  L  J 

where 

A(s)  and  E(s)  are  the  azimuth  and  elevation  mirror  positions,  represented  by  the 
electrical  signals  of  the  position  sensors, 

Vyy(s)  and  Ve(s)  are  voltage  signals  applied  to  the  mirror  drive  system, 

Gaa(s)  and  Gee(s)  are  transfer  functions  describing  the  main  dynamic  channels  of 
the  mirror, 

Gea(s)  and  Gae(s)  are  transfer  functions,  representing  the  cross-effects  between 
the  main  dynamic  channels  of  the  mirror  (so-called  cross-coupling  terms),  and 

argument  s  is  the  Laplace  variable. 

The  following  tasks  were  performed: 

1)  a  5  Hz  square  wave  voltage  signal  with  a  magnitude  of  0.05  V  was  applied  to 
the  azimuth  input  of  the  mirror  drive.  Responses  A(s)=Gaa(s)Va(s)  and 
E(s)=Gea(s)Va(s)  to  this  signal  were  recorded.  The  choice  of  the  5  Hz  frequency  of  the 
input  signal  provided  sufficient  time  for  these  responses  to  reach  steady  state.  The 
responses  contain  complete  information  for  the  estimation  of  transfer  functions  Gaa(s) 
and  Gea(s). 

2)  a  5  Hz  square  wave  voltage  signal  with  a  magnitude  of  0.05  V  was  applied  to 
the  elevation  input  of  the  mirror  drive.  Responses  E(s)=Gee(s)Ve(s)  and 

A(s)  =Gae(s)Ve(s)  to  this  signal  were  recorded.  The  responses  contain  complete 
information  for  the  estimation  of  transfer  functions  Gee(s)  and  Gae(s). 

3)  experiments  1)  and  2)  were  repeated  with  a  step  magnitude  of  0.025  V  to 
detect  possible  nonlinearities  in  the  mirrors  characteristics.  This  is  done  to  insure  that 
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we  are  operating  in  a  linear  portion  of  the  mirrors  response. 

4)  experiments  1),  2),  and  3)  were  conducted  several  times  to  assure  the 
repeatability  of  the  results. 

The  following  figures  represent  the  recorded  responses  of  the  mirror  system. 


Figure  1-1.  Step  response  of  the  azimuth  channel  to  the  azimuth  input. 


Figure  1-2.  Step  response  of  the  elevation  channel  to  the  elevation  input. 


10 


Figure  1-4.  Step  response  of  the  azimuth  channel  to  the  elevation  input. 

Analysis  of  the  obtained  responses  indicates  that 

•  dynamics  of  Az-to-Az  and  El-to-El  channels  is  comprised  of  "main"  dynamics 
and  vibration  (bending  modes).  The  main  dynamics  are  typical  for  a  "one  zero  -  two 
complex  poles  system".  The  bending  modes  have  a  dominant  frequency  of  approximately 
1  KHz  and  settle  within  .015  seconds 

•  in  addition,  the  elevation  response  exhibits  a  100  Hz  mode  with  a  very  low 
magnitude 
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•  the  response  of  the  El-to-Az  and  Az-to-El  channels  exhibits  a  1  KHz  vibration 
with  a  settling  time  of  approximately  .02  sec 

•  magnitudes  of  responses  to  .025  V  signals  are  approximately  50%  lower  than 
magnitudes  of  responses  to  .05  V  signals  which  indicates  that  within  this  range  of 
magnitudes  the  system  is  linear. 

The  following  are  the  transfer  functions  describing  appropriate  dynamic  channels 
of  the  mirror  system  were  established  by  analyzing  of  the  above  responses  of  the  mirror 
systems, 

1026(s+1006)  18s(l+. 000927s) 

^  ^  s2+8000s+l.le7  s2+550s+3.e7(l+000027s) 

1039(s+1100)  20s  8s 

Gpp(s)= - - -  + -  +  - - 

s2+8000s+l.le6  s^+460s+3.e7  s^+375s+8.e5 


100s 


Gae(s)  “ 


s2+260s+3.e7 


Gea(s)— 


120s 


s"+8000s+l.le6 


(1.2) 


These  transfer  functions  represent  the  "main"  dynamics  of  the  Az-to-Az  and  El-to-El 
channels,  and  the  "undesirable"  dynamic  components,  known  as  the  bending  modes. 
Implemented  in  the  Matlab™  software,  transfer  functions  (1.2)  were  used  to  simulate  Az 
and  El  responses  of  the  mirror  position.  Parameters  of  (1.2)  were  manually  tuned  to 
achieve  close  resemblance  of  the  recorded  responses.  Simulated  step  responses  are  given 
in  Figure  1-5  through  Figure  1-8  below. 

The  obtained  results  indicate  that  5  KHz  mirror  with  the  existing  control  system 
has  the  following  flaws: 

•  "unwanted"  high  frequency  oscillatory  transient  component  (bending  modes), 

•  cross  coupling  between  the  Az  and  El  channels  in  the  frequency  range  of  the 
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bending  mode, 

•  nonlinear  frequency  response  within  the  operational  frequency  range, 

•  its  bandwidth  is  definitely  below  the  required  5  Khz. 


Figure  1-5.  Simulated  step  response  of  the  elevation  output  to  the  elevation  input 


•.a 


Figure  1-7.  Simulated  step  response  of  the  elevation  output  to  the  azimuth  input. 


Figure  1-8.  Simulated  step  response  of  the  azimuth  output  to  the  elevation  input. 

B.  Detailed  model  of  a  high  bandwidth  fine  steering  mirror.  The  simple  model  of 
the  high  bandwidth  fine  steering  mirror  was  not  able  to  accurately  predict  the  systems 
response  to  the  addition  of  the  model  reference  controller.  This  led  to  the  development 
of  a  model  with  higher  fidelity  [21].  Although  the  use  of  the  position  sensor  allowed  for 
easier  measurement  of  the  system,  the  optical  quadrant  detector  was  designated  as  the 
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position  sensor.  This  was  justified  by  the  fact  that  in  practice  the  quad  cell  is  the  sensor 
actually  used  to  close  the  position  control  loop. 

The  model  was  developed  based  on  a  series  of  measurements  of  the  output  of  the 
quad  cell  in  response  to  a  step  applied  directly  to  the  mirror  drive  amplifier.  These  were 
an  open  loop  measurements  of  the  actuator  drive  amplifier,  mirror,  and  quad  cell  only. 
The  control  electronics  were  then  modeled  based  on  the  component  values.  The 
following  proceedures  were  followed: 

1)  The  tracking  laser  from  remote  terminal  #2  was  energized. 

2)  The  beam  was  centered  on  the  quadrant  tracking  detector  using  the  coarse 
steering  mirrors  on  the  primary  station  and  by  bias  voltages  applied  to  the  fast  steering 
mirror  drive  drcuitiy 

3)  a  0.25  Hz  square  wave  voltage  signal  with  the  magnitude  of  900  //V  was  applied 

to  the  drive  of  the  Ml/3  actuators.  Recall  that  the  actuators  are  located  in  the  comers 
of  the  miitors  while  the  position  is  measured  relative  to  the  upper-lower  center  of  the 
mirror  (elevation)  and  the  left-right  center  of  the  mirror  (azimuth).  This  corresponds  to 
a  45°  coordinate  rotation.  Responses  A(s)=  ^mwa(®)^mi/3(s)  E(s)— GMiy3E(s)V^jy3(s) 

to  this  signal  were  recorded.  The  choice  of  0.25  Hz  frequency  of  the  input  signal 
provided  sufficient  time  for  these  responses  to  reach  steady-state.  These  responses 
contained  complete  information  for  the  estimation  of  transfer  functions  ^mwaCs)  ^nd 
Gmi/3e(s)- 

4)  a  0.25  Hz  square  wave  voltage  signal  with  the  magnitude  of  900  /zV  was 
applied  to  the  drive  of  the  M2/4  actuators.  Responses  El(s)  “^M2/4e(s)Vm2/4(s) 
Az(s)=G^j2/4a(s)Vm2/4(s)  to  this  signal  were  recorded.  These  responses  contained 
complete  information  for  the  estimation  of  transfer  functions  Gm2/4e(s)  and  Gm2/4a(s)- 

5)  experiments  3)  and  4)  were  repeated  with  the  magnitude  of  1.3  mV  in  order  to 
detect  possible  nonlinearities  in  mirror  characteristics. 

6)  experiments  3),  4),  and  5)  were  conducted  several  times  to  assure  the 
repeatability  of  the  results. 

The  recorded  step  responses  are  given  below.  Figure  1-9  shows  the  response  of 
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the  azimuth  position  to  the  900  //V  step  signal  applied  to  the  Ml/3  driver.  Figure  1-10 
shows  the  response  of  the  elevation  position  to  the  900  /zV  step  signal  applied  to  the 
M2/4  input.  Similar  results  were  obtained  for  the  magnitude  of  the  input  signal  of  1.3 
mV. 


Figure  1-9.  Measured  azimuth  response  to  voltage  applied  to  actuators  Ml/3. 
The  elevation  response  to  the  Ml/3  input  exhibits  a  settling  time  of  about  0.95 
second  with  an  oscillation  frequency  of  24  Hz  with  step  response  that  is  negative  with 
respect  to  the  input. 


time  in  seconds 


Figure  1-10.  Measured  elevation  response  to  voltage  applied  to  actuators  M2/4. 
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The  azimuth  response  exhibits  a  0.6  second  settling  time  and  an  oscillation 
frequency  of  25  Hz.  The  elevation  response  to  the  M2/4  input  exhibits  a  settling  time  of 
about  0.25  second  and  a  oscillation  frequency  of  25  Hz. 

Analysis  of  the  results  once  again  reveals  a  response  that  is  consistent  with  a  one 
zero  and  two  complex  pole  system.  The  same  proceedure  of  hand  tuning  the  transfer 
functions  with  the  MatLab  software  was  used.  The  resultant  transfer  functions  are: 

754(s+7200) 

Gm2/4a(s)  -  "2"'3;3"+ 2465Y 


Gm@/4e(s)  ~ 


-900(s+3200) 

*^+5s'+24655 


-800(s+6907) 
Gmi^a(s)  -  YTYYs + 2465'3 


-900(s+3248) 

^m‘^-(^)=Y4'8sT24665 


(1.3) 


Implemented  in  software,  (1.3)  were  used  to  obtain  the  following  simulated  system 
responses. 


Figure  1-11.  Simulated  response  to  excitation  of  the  M2/4  actuators. 
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Figure  1-12.  Simulated  response  to  excitation  applied  to  the  Ml/3  actuators. 

Once  the  mirror  quad  cell  model  was  complete,  the  control  system  model  was 
established  based  on  the  designed  component  values  and  inherent  nonlinearities.  The 
resulting  overall  model  is  shown  as  a  MATLAB  Simulink  block  diagram  in  Figure  1-13. 
The  designations  for  the  transfer  functions  listed  in  Figure  1-13  relate  them  to  the 
operational  amplifiers  shown  on  the  schematics.  Operational  amplifiers  are  prefixed  with 
a  U  for  example  UlO  and  I  added  the  sufix  letter  to  designate  the  appropriate  channel 
(azimuth  or  elevation).  The  saturation  blocks  were  added  to  model  the  real  saturation 
experienced  by  the  op  amps.  Saturation  occurs  a  little  below  the  op  amp  supply  voltage 
of  ±15V. 

The  system  was  subjected  to  a  ±200  mv  square  wave  with  a  frequency  of  10  Hz. 

A  comparison  of  the  input  to  output  exhibited  a  90  /iS  delay.  The  delay  is  implemented 
in  simulink  through  the  use  of  the  Fade  approximation. 

The  ±200  mv  square  wave  was  introduced  in  the  system  at  the  mirror  controller 
via  the  tracking  offset  port.  The  input  of  this  signal  is  applied  to  U18.  The  output  of 
the  quad  cell  is  measured  at  the  tracking  position  connector  which  is  connected  to  the 
output  of  U19.  Both  connectors  are  BNC  type  and  located  on  the  front  pannel  of  the 
mirror  controller. 
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Figure  1-13.  Simuliuk  representation  of  detailed  model. 


Figure  1-14  and  Figure  1-15  show  the  systems  measured  response  to  a  square 
wave  applied  to  the  elevation  and  azimuth  input.  The  data  was  collected  using  a  Phillips 
digital  oscilloscope  and  the  trace  includes  a  500  /iS  delay.  Both  channels  exhibit  a  25% 
overshoot  in  the  principal  channel  and  some  cross  channel  coupling.  Notice  that  this 
coupling  is  smaller  than  presented  for  the  case  where  the  DIT  position  sensors  were  used 
to  close  the  position  control  loop. 


Figure  1-16  and  Figure  1-17  show  the  simulated  response  to  a  step  input  for  both 
the  elevation  and  azimuth  channel. 


Figure  1-16.  Simulated  response  to  a  step  applied  to  the  elevation  channel. 


Figure  1-17.  Simulated  response  to  a  step  applied  to  the  azimuith  channel. 

12  Piezo-electric  mirror 

Piezo-electric  mirror  intended  for  fine  steering  of  the  laser  beam  was  a  Physik 
Instrumente  S-330  Low  voltage,  3  axis,  tilting  mirror.  The  mathematical  model  of  this 
mirror,  including  its  control  circuitiy,  is  based  on  the  series  of  transient  responses 
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recorded  during  laboratory  testing.  A  square  wave  is  applied  to  each  of  the  input 
channels  and  the  output  of  the  quadrant  detector  is  recorded.  A  mathematical  model  of 
the  system  is  described  as  follows: 


r 

r 

1  r  1 

A{S) 

|Gai(s) 

=  1 

(s) 

Ga3(s)  I  |v,(s)  1 

1  iv^Cs)  1 

E(S) 

|Gbi(s) 

G„(s) 

G«(S)  1  |V3(S)  1 

L 

J  L  J 

Once  again  the  following  tasks  were  performed: 

1)  a  1  Hz  square  wave  voltage  signal  with  a  magnitude  of  0.02  V  was  applied  to 
the  channel  1  input  of  the  mirror  drive.  Responses  A(s)=Gai(s)Vi(s)  and 
E(s)=Gei(s)Vi(s)  to  this  signal  were  recorded.  The  choice  of  the  5  Hz  frequency  of  the 
input  signal  provided  sufficient  time  for  these  responses  to  reach  steady  state.  The 
responses  contain  complete  information  for  the  estimation  of  transfer  functions  Gai(s) 
and  Gei(s). 

2)  a  1  Hz  square  wave  voltage  signal  with  a  magnitude  of  0.02  V  was  applied  to 
the  channel  2  input  of  the  mirror  drive.  Responses  A(s)=Ga2(s)V2(s)  and 
E(s)=Ge2(s)V2(s)  to  this  signal  were  recorded.  The  responses  contain  complete 
information  for  the  estimation  of  transfer  functions  Gy^(s)  and 

3)  a  1  Hz  square  wave  voltage  signal  with  a  magnitude  of  0.02  V  was  applied  to 
the  channel  3  input  of  the  mirror  drive.  Responses  A(s)=Ga3(s)V3(s)  and 
E(s)=Ge3(s)V3(s)  to  this  signal  were  recorded.  The  responses  contain  complete 
information  for  the  estimation  of  transfer  fimctions  Ga3(s)  and  Ge3(s). 

4)  experiments  1)  ,  2)  and  3)  were  repeated  with  a  step  magnitude  of  0.04  V  to 
detect  possible  nonlinearities  in  the  mirrors  characteristics.  This  is  done  to  insure  that 
we  are  operating  in  a  linear  portion  of  the  mirrors  response. 

5)  experiments  1),  2),  and  3)  were  conducted  several  times  to  assure  the 
repeatability  of  the  results. 

Figures  1-18,  1-19,  1-20  below  represent  the  measured  step  responses  of  channels 
1,  2,  and  3  respectivly.  These  measurements  reveal  a  fairly  significant  degree  of  cross 
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coupling,  or  exchange  of  bending  moodes  as  evidenced  by  the  beat  present  in  the 
channel  2  and  3  response.  These  responses  are  consistent  with  a  3  zero,  4  pole  response. 


Figure  1-18.  Step  response  for  channel  1. 


Figure  1-19.  Step  response  for  channel  2 


-21- 


Figure  1-20.  Step  response  channel  3. 

A  simulation  model  was  developed  on  the  basis  of  measurements  and  then  hand 
tuned  until  it  resembled  the  above  responses.  Figures  1-21,  1-22,  and  1-23  depict  the 
simulated  responses  of  channels  1,2  and  3  of  the  model.  Figure  1-24  represents  resultant 
model  in  the  form  of  a  simulation  setup. 


Figure  1-21.  Simulated  channel  1  step  response. 


Figure  1-22.  Simulated  channel  2  step  response. 


Figure  1-23.  Simulated  channel  3  step  response. 
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Figure  1-24.  Simulink  model  of  piezo  mirror. 


-24- 


2.  MODEL  REFERENCE  CONTROL  OF  BEAM  STEERING  MIRRORS 

A  mechanical  beam  steering  element  is  the  critical  component  of  the  entire  PAT 
system,  responsible  for  its  major  operational  limitations.  Much  work  has  been  devoted 
to  the  development  of  high  performance/bandwidth  steering  mirrors  and  other  non¬ 
mechanical  beam  steering  devices.  Currently,  these  devices  are  available  only  as 
developmental  models  and  consequently  are  veiy  expensive.  Commonly,  a  high 
bandwidth  feedback  control  is  used  to  stabilize  mechanical  beam  steerers. 

Applications  of  advanced  control  techniques,  such  as  model  reference  (MR) 
control,  applied  to  commercially  available  mirrors  provides  a  means  of  achieving  the 
required  performance  at  a  reasonable  cost.  MR  approach  implies  that  the  response  of 
the  existing  physical  system  is  "forced"  to  follow  the  response  of  a  model  implemented  in 
analog  circuitry  or  in  software.  The  model  is  subjected  to  the  same  commanded  input  as 
the  laser  beam  positioning  system,  but  the  model  response  is  "perfect"  and  not  affected 
by  satellite  jitter  or  any  other  undesirable  effects.  MR  approach  was  suggested  for 
extending  performance  limits  of  an  optical  mirror  and  improving  its  dynamics  [12,13,16]. 
While  the  MR  scheme  eliminates  any  discrepancy  between  the  responses  of  the  physical 
system  and  the  model,  it  has  the  potential  for  jitter  rejection.  As  with  any  physical 
system  subjected  to  long  term  effects  of  adverse  operational  conditions  and  natural  wear, 
performance  of  a  beam  positioning  system  is  expected  to  deteriorate  with  time.  The 
required  system  dynamics,  represented  by  the  simulation  code  of  a  reference  model,  is 
time-invariant.  Therefore,  the  required  performance  characteristics  of  a  beam 
positioning  system  can  be  assured  and  maintained  via  application  of  MR  approach.  It  is 
important  that  MR  technology  does  not  require  any  additional  hardware  and  numerically 
is  much  simpler  than  the  feedforward  adaptive  jitter  compensation. 

MR  approach  is  suggested  for  further  improvement  of  the  existing  pointing- 
acquisition-tracking  (PAT)  systems  intended  for  laser  communication:  reduction  of 
adverse  effect  of  satellite  jitter,  and  increase  of  system  robustness.  It  can  be 
implemented  via  a  number  of  schemes.  Inherent  limitations,  advantages  and 
disadvantages  of  the  particular  schemes,  and  various  aspects  of  practical  implementation 
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of  this  approach  are  investigated  below. 

A  typical  beam  steerer  includes  a  high  performance  optical  mirror,  drives,  position 
sensors,  and  feedback  control  circuitry.  The  mirror  has  two  control  channels,  elevation 
and  azimuth.  The  beam  position,  defined  by  a  quadrant  detector,  is  controlled  by  the 
application  of  the  appropriate  control  efforts  in  compliance  with  the  reference  signals. 
While  it  is  desired  that  dynamic  channels  be  controlled  independently,  the  cross-coupling 
effects,  which  can  be  viewed  as  the  "exchange  of  bending  modes"  could  be  always 
observed.  Analyses  of  closed-loop  dynamics  of  particular  channels  indicate  nonlinearity 
and  quite  limited  bandwidth  of  the  appropriate  frequency  responses.  In  addition, 
satellite  jitter  constitutes  a  permanent  source  of  beam  positioning  errors. 

Application  of  MR  control  schemes  has  the  potential  for  alleviation  of  the  above 
listed  problems.  Conceptually,  it  implies  that  a  reference  model  representing  the 
desired  closed-loop  system  performance  be  established  in  the  form  of  an  analog  circuit 
or  a  computer  code.  In  a  MR  scheme,  an  additional  control  signal  defined  by  a  specially 
designed  controller,  "forces"  the  beam  steerer  to  follow  the  response  of  the  reference 
model.  Understandably,  a  reference  model  not  only  has  perfect  dynamic  characteristics, 
but  cannot  be  affected  by  satellite  jitter  and  hardware  deterioration.  Therefore,  in 
addition  to  improved  dynamics,  this  approach  results  in  jitter  rejection  and  increased 
system  robustness.  Limitations  of  the  MR  approach  are  related  to  the  limited  magnitude 
of  control  effort  which  can  be  applied  to  a  beam  steerer,  and  the  complexity  of  its 
implementation. 

Conventional  MR  schemes,  capable  of  decoupling,  extended  system  bandwidth, 
reduction  of  jitter  effects,  and  high  system  robustness,  require  state-variable 
implementation  and  are  overly  complex  (numerically  intensive)  for  many  applications.  In 
many  instances,  when  some  of  the  above  listed  effects  are  not  required,  particular 
simplifications  of  a  "general"  MR  control  can  be  suggested. 

Described  below,  is  a  model-following,  rather  than  MR,  control  scheme  utilizes 
so-called  System  Error  Equation  Method.  This  method  results  in  performance 
modification  and  jitter  rejection,  but  not  in  robust  performance.  Such  an  approach  can 
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be  recommended  when  system  dynamics  is  fairly  well  defined  and  assumed  to  be  time- 
invariant.  The  system  utilizes  readily  available  input-output  signals. 

The  second  system,  described  herein,  utilizes  an  adaptive  predictor  and  is 
intended  for  digital  implementation.  It  exhibits  all  features  of  a  MR  scheme,  but  does 
not  require  a  state  observer. 

2.1  System  Error  Equation  Method  of  Model-Following 

Consider  closed-loop  dynamics  of  a  mechanical  beam  steerer  described  in  the 


form. 


Y(s)=G(s)Up(s) 

the  reference  model  is  defined  as  follows. 

(2.1) 

Ym(s)=M(s)U(s) 

the  output  error  is  defined  as. 

(2.2) 

E(s)=Ym(s)-Y(s) 

and  the  controller  implementing  the  MR  concept  is. 

(2.3) 

Up(s)=U(s)-hW(s)E(s) 

(2.4) 

where 


G(s)=[gij(s),  i,j=l,2]  and  M(s)=[mij(s),  i,j=l,2]  are  appropriate  transfer  matrices; 
it  is  expected  that  G(s)  is  non-diagonal,  and  M(s)  is  a  diagonal  matrix, 

Y=[yi  yj]’,  and  YM=[yMi  yM2]’  ^^e  elevation  and  azimuth  position  components  of 
the  existing  beam  steerer  and  the  reference  model, 

Up=[upi  Upj]’,  and  U=[ui  U2]’  are  control  signals  applied  to  the  existing  beam 
steerer  and  to  the  reference  model  (also  known  as  the  commanded  input),  and 
W(s)=[Wij(s),  i,j=l,2]  is  the  MR  controller. 

Rewrite  equation  (2.1)  as  follows, 

Y(s)=G(s)U(s)-hG(s)W(s)E(s)=G(s)U(s)-hG(s)W(s)[YM(s)-Yp(s)]= 

G(s)U(s)-fG(s)W(s)Ym(s)-G(s)W(s)Yp(s), 

then 
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Y(s)+G(s)W(s)Yp(s)=G(s)U(s)+G(s)W(s)Ym(s),  or 
[I + G(s)  W(s)]  Y(s) = [G(s) + G(s)W(s)M(s)]U(s) 

Finally, 

Y(s)  =  [I + G(s)W(s)]'G(s)[I + W(s)M(s)]U(s)  (2.5) 

Redefine  the  output  error  as 

E(s)={M(s)  -  [I+G(s)W(s)]  ‘G(s)[I+W(s)M(s)]}U(s)  (2.6) 

Define  an  arbitrary  transfer  matrix 

Q(s)=[qij(s),  i,j=l,2]  (2.7) 

representing  the  desired  dynamics  of  the  "commanded  input  -  error"  channels.  The 
implication  of  the  word  "desired"  are, 

-  transfer  matrix  Q(s)  must  be  stable, 

-  settling  time  not  to  exceed  1/2  of  the  settling  time  of  the  reference  model,  and 

-  sufficiently  small  steady-state  values  of  the  error  response,  which  can  be  assured 
by  l|Q(0)||<  1 

Note,  that  there  is  no  need  to  achieve  the  error  decoupling.  Then  the  following  equation 
should  be  established, 

M(s)  -  [I-l-G(s)W(s)]-‘G(s)[H-W(s)M(s)]=Q(s)  (2.8) 

which  results  in  the  definition  of  the  controller  W(s), 

W(s)=G(s)-iM(s)Q(s)  '  -  G(s)-*  -  Q(s)-‘  (2.9) 

It  should  be  noted  that  the  major  advantage  of  this  approach  is  a  veiy  high  flexibility  in 
the  particular  choice  of  error  equation  matrix  Q(s).  Indeed,  upon  exploration  of 
equation  (9)  allows  for  the  choice  of  such  Q(s)  which  would  result  in  the  sufficiently 
simple  definition  of  the  controller  W(s).  The  disadvantage  of  this  approach  is  quite 
typical  of  any  linear  model-following  design:  it  implies  that  matrix  G(s)  is  accurately 
known  and  time-invariant.  However,  according  to  our  preliminaiy  results,  the  developed 
model-following  scheme  can  be  successfully  extended  to  a  MR  scheme  utilizing  input- 
output  signals. 

The  following  numerical  example  is  based  on  a  simplified  mathematical 
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description  of  a  beam  steerer.  Particular  elements  of  the  matrix  G(s)  are, 

gii(s)=2.6186E8/(s'+444s+2.1822E8) 

gi2(s)=2.0731E8/(s2+444s+2.1822E8) 

g2i(s)  =  l.E8[1.005s^+428.39s+ 1.3735E8]/[(s^+471s+ 1.5824E8)(s^+400s+ 1.24E8)] 
g22(s)=-l.E8[1.067s2+463.16s+1.4985E8]/[(s2+471s+1.5824E8)(s^+400s+1.24E8)]  (2.10) 
The  reference  model  transfer  matrix  is  defined  as, 
mii(s) = m22(s) = 3000/(s + 3000) 

mi2(s)=m2i(s)=0 

The  following  is  the  definition  of  the  "error  transfer  matrix"  Q(s), 
qii  (s) = q22(s) = .001  *s/(s + 6000) 

qi2(s)=q2i(s)=0 

leading  to  the  definition  of  controller  W(s), 

Wii(s)=-l.EllWo(s)[2.6s*+1.925E4s'+1.87E8s2+1.23E12+1.49E15] 

Wi2(s)=37.1wo(s)[(s^+426s+1.37E8)(s^+444s+2.18E8)(s2+400s+1.14E8)] 

W2i(s)=76.5wo(s)[(s2+471s+1.58E8)(s2+400s+1.14E8)(s2+400s+1.24E8)] 

W22(s)=-1.E11wo(s)[3.24s'+3.29E4s’+2.72E8s'+1.86E12+5.09E15] 

where 

Wo(s)=[s(s+3000)(s*+831s^+2.56E8s"+1.059Ells+1.62E16]-' 

Computer  simulation  results,  see  pages  30,  31,  indicate  that  the  obtained 
controller,  even  after  reasonable  simplification,  is  capable  of  "enforcing"  the  required 
i  system  dynamics  on  the  beam  steerer  and  reduced  the  effect  of  satellite  jitter  on  the 

beam  positioning  accuracy. 
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STEP  RESPONSE  OF  THE  ORIGINAL  S VS TEN  TO  THE  AZINUTH  INPUT 


STEP  RESPONSE  OF  THE  ORIGINAL  SYSTEM  TO  THE  ELEUATION  INPUT 
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12  System  Error  Equation  Method.  State-Variable  Formulation 

Consider  a  controlled  plant  and  a  reference  model  and  defined  by  state  equations: 

X’=ApX+BpUp  (2.11) 

Z’=AmZ+BmU  (2.12) 

where 

X=[xi,  X2,...,x„f  and  Z=[zi,  Z2,.  -,Znr  transposed  state  vectors  of  the  controlled 
plant  and  the  model, 

U  and  Up  are  vectors  of  reference  signals  (commanded  input)  and  control  effort 
applied  to  the  plant  and  the  model  through  matrices  Bp  and  B,^j, 

Ap  and  A^  are  fundamental  matrices  of  the  controlled  plant  and  the  reference 
model,  and 

^  is  transpose  symbol. 

Define  the  state  generalized  error  as: 

E=Z-X  (2.13) 

The  control  law  is  defined  as  follows: 

Up=-KpX+KuU-l-WE  (2.14) 

where 

Kp  represents  the  existing  state-variable  feedback  controller, 

Ku  is  the  matrix  gain  in  the  reference  channel,  and 

W  is  the  matrix  representing  the  model  reference  controller,  responsible  for  the 
elimination  of  the  discrepancy  between  the  state  vectors  of  the  plant  and  the  model. 
Redefine  the  plant  with  its  state-variable  feedback: 

X’=ApX-BpKpX-l-BpKuU-l-BpWE  (2.15) 

Obtain  and  transform  the  expression  of  the  error  derivative: 

E’=Z’-X’,  or 

E’=AmZ-I-BmU-ApX+BpKpX-BpKuU-BpWE,  or 

E’  =  AmZ-I-  BmU-Ap(Z-E)  -I-  BpKp(Z-E)-BpKuU-BpWE,  or 
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E’  =  AmZ+BmU-ApZ+ ApE+BpKpZ-BpKpE-BpKuU-BpWE,  or 

E’  =  ApE-BpKpE-BpWE+ AmZ-ApZ+BpKpZ+BmU-BpKuU,  or 

E’  =  (Ap-BpKp-BpW)E+(AM-Ap+BpKp)Z+(BM-BpKu)U  (2.16) 

Note  that  matrix-vector  equation  (2.16)  describes  dynamics  of  the  state 
generalized  error  driven  by  two  "external"  signals,  commanded  input  U,  and  state  vector 
of  the  reference  model,  Z.  One  should  realize  that  both  signals,  U(t)  and  Z(t),  are  "well 
behaving"  signals,  and  therefore,  the  desired  convergence  rate  of  error  E  is  dependent 
upon  matrix  Q=Ap-BpKp-BpW. 

Introduce  matrix  a  nxn  matrix  Q  such  that, 

-  its  eigenvalues  assure  the  necessaiy  error  convergence  rate, 

•  ||Q''(Aj^-Ap-f  Bpl^)||<  1  and  ||Q  ^(B^-BpKu)  ||<  1 
to  assure  small  steady-state  error  values  of  the  "natural  motion"  terms  of  the  error. 

Then,  the  following  expression  for  the  model  reference  controller  can  be  obtained: 
W=Bp-i(Ap-Q)  -  Kp  (2.17) 

It  can  be  seen  that  the  existence  of  a  solution  requires  that  matrix  Bp  be  non-singular. 
Numerical  example. 

Assume  the  following  definition  of  the  controlled  plant  and  the  reference  model: 
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One  can  realize  that  the  controlled  plant  is  coupled,  and  the  reference  model  is 
uncoupled  and  reflects  the  following  design  specifications:  overshoot  of  0%  and  settling 
time  of  .5  sec.  The  following  state-variable  controller  was  introduced  in  the  system: 

r  •] 

I  248  105  21  0  0  0| 

K.=  0  0  0  251  107  20  I 

L  J 

which  indicates  that  the  designers  of  the  original  system  were  interested  in  the  overshoot 
of  0%  and  settling  time  of  1  sec  for  each  dynamic  channel.  The  matrix  gain  in  the 
reference  channel  of  this  system  was  defined  as 

r 

250  0 

Ko=  0  250 

Simulated  step  responses  of  the  designed  system  are  shown  on  page  35.  They  indicate 
that  the  designed  system  exhibits  cross-coupling  between  channels.  However,  in  the 
situation  when  parameters  of  the  plant  are  not  accurately  known,  it  could  be  expected 
that  system  performance  will  be  much  more  different  from  the  original  design 
specifications. 

Introduction  of  the  MR  controller  is  expected  to  significantly  modify  system 
performance  and  reduce  the  cross-coupling  effects.  Following  the  equation  error 
approach,  first  introduce  matrix  Q: 
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which  reflects  the  requirement  that  state  generalized  error  converge  within  .4  sec 

Now  the  MR  controller  represented  by  matrix  W  can  be  defined  according  to 
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equation  (2.17).  First,  realize  that  matrix  Bp  is  not  a  square  one  and  its  inverse  in 
(2.17)  must  be  replaced  by  a  left  pseudo-inverse  H=(Bp'^Bp)'^Bp'r.  Note  that 

HBp = (BpTBp)-'BpTBp = (BpTBp)-'  (Bp^Bp) = I 
providing  that  the  inverse  (Bp^Bp)'*  exists.  In  our  situation, 

r  1 

I  .0  .0  .609756  .0  .0  .487805| 

H=l  .0  .0  -.487805  .0  .0  .609756| 

L  J 


and,  finally  matrix  W  is  defined  as: 

r  1 

I  9.5678e3  632.805  17.4146  7. 854163  591.219  30.2439 
W=| -7. 852763  -590.244  -30.7317  9.566763  632.024  17.8049 

L 

Simulated  step  responses  of  the  designed  MR  system  are  shown  on  page  37.  They 
indicate  that  the  designed  system  exhibits  modified  performance  and  cross-coupling 
transients  decay  within  the  time  designated  for  the  error  convergence. 


2J  Decoupling  MR  Control  with  an  Adaptive  Predictor 
Consider  a  beam  steerer  described  by  a  transfer  matrix  G(s).  Assume  that 
elements  of  the  same  row  of  matrix  G(s)  are  brought  to  a  common  denominator. 

Convert  matrix  G(s)  to  its  discrete-time  equivalent,  G(z),  assuming  that  the  steerer  is 
driven  through  a  zero-order-hold.  Then 

the  following  difference  equations  could  be  used  for  the  description  of  the  beam  steerer, 
yi(i)=-Saijyi(i-j)  -I-  SbijUpi(i-j)  +  ScijUp2(i-j) 

y2(i)=-Sa2jyi(i-j)  +  Sb2jUpi(i-j)  +  Sc2jUp2(i-j),  (2.18) 

where  i  is  the  discrete-time  index.  One  can  see  that  summation  in  (11)  is  performed 
with  respect  to  time  index  j= 1,2,3,...  ,  and  coefficients  a^j,  bjj,  c^,  a2j,  bjj,  and  Cjj,  are  the 
parameters  of  the  z-domain  transfer  functions  gii(z),  £12(2),  £21(2).  and  £22(2),  and  as  such 
are  dependent  on  the  time  discretization  step. 

Note  that  since  y,, ,  k=l,2  are  azimuth  and  elevation  positions  of  the  laser  beam 
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accurately  defined  by  the  quadrant  detector,  and  Up^ ,  k=l,2,  are  known  control  signals, 
parameters  ajj,  bjj,  Cjj,  a2j,  b2j,  and  Cjj,  can  be  continuously  updated  by  the  application  of  a 
finite -memoiy  recursive  least  squares  procedure  [20].  Although  the  success  of  the 
parameter  updating  can  be  cissured  by  good  initial  values  of  these  parameters  and  low 
measurement  noise,  it  is  known  that  when  required,  the  procedure  can  be  enhanced  by 
orthogonalization  of  variables  yi(i-j),  y2(i-j),  Upi(t-j)»  and  Up2(i-j)  and  consequent 
modification  of  the  obtained  estimates.  Implementation  of  the  parameter  updating 
procedure  is  instrumental  in  situations  when  parameters  of  the  controlled  process  are 
poorly  known  and/or  time-dependent  due  to  hardware  deterioration  and  environmental 
effects. 

Availability  of  accurately  estimated  parameters  of  the  equations  (2.18)  facilitates 
the  use  of  these  equations  for  the  prediction  of  azimuth  and  elevation  positions.  Then 
continuous  parameter  updating  provides  the  assurance  of  the  prediction  accuracy. 

Indeed,  if  (2,18)  holds  then  the  following  "one  step  ahead"  prediction  can  be  achieved  as 
yi(i-l-l)=-Saijyi(i-j+l)  +  SbijUpi(i-j+l)  -I-  SCijUp2(i-j+l) 

y2(i+l)=-Sa2jyi(i-j+l)  +  Sb2jUpi(i-j-f-l)  -I-  Sc2jUp2(i-j+l)  (2.19) 

Introduce  a  reference  model  in  the  form  of  two  fixed-parameter  difference  equations, 
reflecting  the  desired  closed-loop  dynamics  of  the  beam  steerer, 
yMi(i+l)=-SmijyMi(i-j)  +  SnijUi(i-j) 

yM2(i+l)=-Sm2jyMi(i-j)  +  Sn2jUi(i-j)  (2.20) 

One  can  see  that  equations  (2.20)  represent  a  decoupled  system.  Computation  of  the 
predicted  azimuth  and  elevation,  yi(i+l)  and  y2(i-l-l),  and  their  reference  values 
yMi(i+l)*  yMzO+l)  facilitates  computation  of  the  expected  position  errors, 

ei(i-f-l)=yMi(i+l)-yi(i+l)  (2.21) 

e2(i+l)=yM2(i+l)-y2(i+l) 

Define  increments  of  control  signals,  AUpi(i)  and  AUp2(i)  to  assure  that  position  errors 
have  zero  values, 

biiAupi(i)  +  qiAup2(i)  =  -ei(i-»-l) 
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(2.22) 


b2iAupi(i)  +  C2iAup2(i)  =  -ejCi+l) 

Then  control  signals  Upi(i+1)  and  Up2(i+1)  can  be  defined  as  follows, 

Upi(i + 1) =Upi(i + 1)  ^^pi(i) 

Up2(i+l)=Up2(i+l)+^Up2(i)  (2.23) 

The  proposed  MR  control  procedure  performing  parameter  updating,  solution  of 
equations  (2.22),  and  definition  of  the  control  signals  as  per  (2.23)  at  each  step  of  the 
discrete-time  procedure,  i= 1,2,3,...,  has  been  implemented  in  software  and  tested  by 
simulation.  It  is  important  that  the  proposed  decoupling  MR  technique  can  be  extended 
from  two  to  any  practical  number  of  input-output  channels. 

The  following  are  difference  equations  representing  the  beam  steerer,  originally 
defined  by  transfer  functions,  at  sampling  frequency  32.5  KHz, 

yi(i) = 1.7857yi(i-l)-.9865yi(i-2) + .1207[upi(i-l) + .9954upi(i-2)]+ 

.0956[up2(i-l) + .9954up2(i-2)] 

y2(i)=3.7122y2(i-l)-5.41815y2(i-2)+3.6628y2(i-3)-.9736y2(i-4)-t- 

.0466[upi(i-l)-.8649upi(i-2)-.8649upi(i-3)  -I-  .9825upi(i-4)]  - 

.0495[up2(i-l)-.8612up2(i-2)-.8616up2(i-3)-l-.9823up2(i-4)]] 

The  reference  model  is  represented  by  expressions 
y^i(i)=.912yMi(i-l)  +  ,08799ui(i-l) 
y^2(i)=-912yM2(i-l)  +  .08799u2(i-l) 

Computer  simulation  results  indicate  that  recursive  least  squares  procedure  allows 
for  fast  and  reliable  convergence  of  the  parameter  updating  process,  both  in  the 
estimation  and  tracking  modes.  The  control  procedure,  especially  upon  convergence  of 
the  parameter  updating  process,  allows  for  successful  model  following.  Utilization  of  the 
continuous  parameter  updating  provides  a  high  degree  of  robustness  of  the  resultant 

control  system. 

Graphs  of  page  40  indicate  that  a  high  degree  of  decoupling  and  model-following 
has  been  achieved  by  the  application  of  the  above  approach.  A  practical  question, 
however,  arises:  how  crucial  the  accuracy  of  prediction  is  for  the  performance  of  the  MR 
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system?  The  answer  to  this  question  is  provided  by  the  graphs  of  pp.  41  and  42.  During 
the  simulation  experiment,  the  "true"  output  of  the  process  has  been  contaminated  by 
noise  which  is  not  accounted  for  by  the  predictor.  In  spite  that  the  level  of  noise  is 
approximately  30%  of  the  signal  magnitude,  the  proposed  procedure  allows  for  successful 

decoupling  and  model  following. 

Parameter  tracking.  While  the  necessary  accuracy  of  model-following  is  assured  by 
the  existing  feedback  mechanisms,  the  proposed  control  law  is  dependent  only  on  the 
ability  to  provide  accurate  on-line  estimation/tracking  of  the  parameters  a,  b,  c  of 
equation  (2.19).  This  task  can  be  carried  out  by  the  following  procedure  implementing 
the  Recursive  Least  Squares  Method  (RLSM)  with  exponential  forgetting  [Astrom,  K.J. 
and  Wittenmark,  B.,  "Adaptive  Control,"  Addison-Wesley  Publishing  Company].  In 
order  to  describe  RLSM  application  to  our  problem,  define  first  equation  of  (2.19)  in 
more  detail: 


N  M  K 

yi(i-l-l)=-S  aijyi(i-j+l)  +  2  bijUpi(i-j+l)  +  S  CijUi,20*j‘'‘l)  (2.26) 

j=l  j=l  j=l 

where 

N,  M,  and  K  are  defined  by  the  order  of  denominator  and  numerators  of  the 
appropriate  z-domain  transfer  functions. 


Define  the  following  vectors, 

X(i)=[yi(i-j+l)d=lv..,N,  Upi(i-j-l-l),j=l,...,M,  Up2(i-j-f-l),j=l . K  f 

A=[aij,  j=l,...,N,  bjj,  j=l,...,M,  Cjj,  j— 1,...,K] 


(2.27) 

(2.28) 


then  the  RLSM  equations  can  be  defined  as  follows. 


K(i-Hl)=P(i)X(H-l)[X-»-X(i-l-l)^P(i)X(i-J-l)]'  (2.29) 

a(i-H)=a(i)+K(i-Hl)[y(i+l)-a(i)’"X(i-H)]  (2.30) 

P(i +1)= [I-K(i)X(i  -I-  l)’"]P(i)/X  (2.31) 


where 


i= 1,2,3,...  discrete  time  index, 

a(i),  i= 1,2,3,...  is  a  series  of  estimates  of  vector  A,  generated  by  the  RLSM 
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procedure, 

vector  X(i+1)  represents  recorded  values  of  the  appropriate  variables  as  per 

(2.27), 

vector  a(i)  represents  the  most  recent  vector  of  parameter  estimates, 

P(i)  and  K(i+1)  are  LxL  matrices,  L=N+M+K  (note  that  square  brackets  of 
(2.29)  contain  a  scaler  and  this  expression  does  not  imply  matrix  inversion), 

y(i+l)=yi(i+l)  is  the  most  recent  value  of  this  variable  recorded, 

A  is  the  forgetting  parameter,  0<  A  <1. 

Successful  implementation  of  a  RLSM  procedure  is  dependent  upon  the  choice  of 
its  initial  parameters  A(0)  and  P(0),  measurement  accuracy,  and  the  value  of  the 
forgetting  factor  A.. 

The  choice  of  A(0)  should  be  consistent  with  the  available  parameter  values  of 
the  appropriate  z-domain  transfer  functions  of  the  system. 

It  is  recommended  that  P(0)  be  defined  as  a  diagonal  matrix,  and  the  value  of.  .1 
should  be  assigned  to  all  diagonal  elements.  This  choice  may  be  reconsidered  depending 
on  the  results  of  the  trial  runs  of  the  RLSM  procedure:  noticeable  oscillations  of  the 
parameter  estimates,  converging  or  diverging,  indicate  that  the  values  of  the  diagonal 
terms  of  matrix  P(0)  should  be  reduced.  Slow  convergence  of  the  parameter  estimates 
indicates  that  the  values  of  the  diagonal  terms  of  matrix  P(0)  should  be  increased. 

The  choice  of  parameter  A  controls  the  balance  between  the  ability  of  the  RLSM 
procedure  to  operate  under  measurement  noise  and  the  ability  to  track  time-vaiying 
parameters  A.  Parameter  A  defines  the  "memoiy  size"  of  the  RLSM  procedure,  i.e. 
determines  the  number  of  most  recent  measurements  which  are  primarily  responsible  for 
the  values  of  estimates  A.  The  memoiy  of  the  RLSM  procedure  is  defined  as  1/(1-A). 
One  can  see  that,  while  0<A<1,  an  increase  of  the  parameter  A  results  in  the  reduction 
of  memory  size  and,  therefore,  improves  tracking  capability  of  the  procedure.  A 
decrease  of  the  parameter  A  results  in  the  increase  of  memoiy  size  and,  therefore, 
improves  the  ability  of  the  procedure  to  "average  out"  measurement  noise.  At  the  same 
time,  tracking  capability  of  the  procedure  decreases. 
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A  FORTRAN  text  of  an  RLSM  procedure  developed  for  this  project  can  be  seen 


below. 


SUBROUTINE  RLSM (N, TH, Y , FAI , P, FF) 
DIMENSION  TH(IO) ,G(10) , FAI (10) ,P(10,10) 
DIMENSION  0(10,10) ,D(10,10) ,H(10,10) 

DO  6  1=1, N 
6  H(1,I)=FAI(I) 

CALL  MULT(H,P,1,N,N,C) 

CALL  MULT (C, FAI, 1,N,1,D) 
E=1.0/(D(1,1)+FF) 

DO  1  1=1, N 

1  Y=Y-FAI  (I) 'THd) 

CALL  MULT(P,FAI,N,N,1,C) 

DO  2  1=1, N 

2  G(I)=C(I,1)*E 
DO  3  1=1, N 

3  TH(I)=TH(I)+G(I) *Y 
DO  4  1=1, N 

DO  4  J=1,N 
D(I,J)=-G(I)*FAI(J) 

IF  (I.EQ.J)  D(I, J)=D(I, J)+1.0 

4  CONTINUE 

CALL  MULT(D,P,N,N,N,C) 

DO  5  1=1, N 
DO  5  J=1,N 

5  P(I,J)=C(I,J)/FF 
RETURN 

END 

C 

SUBROUTINE  MULT(A,B,M,K,N,C) 

C  A(MXK)  B(KXN)  C(MxN) 

DIMENSION  A(10,10) ,B(10,10) ,C(10,10) 

DO  1  1=1, M 
DO  2  J=1,N 
C(I,J)=0.0 
DO  3  KK=1,K 

3  C(I,J)=C(I,J)+A(I,KK)*B(KK,J) 

2  CONTINUE 

1  CONTINUE 

RETURN 
END 


Graphs  given  on  page  46  illustrate  the  ability  of  the  procedure  to  estimate  and  to 
track  system  parameters.  To  dramatize  the  situation,  initial  parameter  values  are  chosen 
to  be  zeros. 

Graphs  of  page  47  illustrate  the  ability  of  the  RLSM  estimates  to  follow  time- 
vaiying  "true"  parameters  of  the  system. 
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Accuracy  analysis  of  the  MR  system.  The  purpose  of  this  section  is  to  assess  the 
risk  of  the  application  of  the  proposed  above  procedure  associated  Avith  the  fact  that 
RLSM  estimates  of  system  parameters  may  or  may  not  become  effectively  close  to  their 
"true"  values.  In  order  to  perform  this  analysis,  rewrite  equations  (2.19)  as  follows 
N  M  K 

yi(i+l)=-Zl  aijyi(i-j+l)  +  S  bijUpi(i-j+l)  +  S  CijUp2(i-j+l)  +biiUpi(i)  +  CiiUp2(i) 

j=l  j=2  j=2 

N  M  K 

y2(i+l)=-S  a2jy2(i-j+l)  +  S  b2jUpi(i-j+l)  +  S  C2jUp2(i-j+l)  +b22Upi(i)  +  C22Up2(i)  (2.32) 

j=l  j=2  j=2 

Introduce  the  following  vectors, 

V(i)=[yi(i-j+l),  y2(i‘j"l'l)>  %i(^‘j"*‘l)»  j~2,...,M,  Up2(i-j+l),  j— 2,...,K] 

U(i)=[upi(i)  Up2(i)f, 

Y(i+l)=[y,(i+l)  y2(i+l)r,  and 
YM(i+l)=[yMi(i+l)yM2(i+l)f 

and  matrices  of  appropriate  dimensions,  R  and  Q,  such  that  equations  (2.32)  could  be 
rewritten  in  the  matrix-vector  form 

Y(i-l-l)=RV(i)+QU(i)  (2.33) 

One  can  see  that  matrix  R  contains  most  of  parameters  of  the  z-domain  transfer  function 
of  the  controlled  plant,  and  matrix  Q  id  defined  as 

^11  I 

Q  =  bja  Cja  I 

L  J 

Assume  that  matrix  <I>(i)  and  matrix  0(i)  are  obtained  from  R  and  Q  correspondingly  by 
replacing  appropriate  ^stem  parameters  by  their  RLSM  estimates.  Then  the  predicted 
system  outputs  can  be  expressed  as 

YPRED(i  ^  <I)(i)  V(i) + 0(i)U(i) 

and  the  predicted  system  error  as 

EPRED(i+l)=YM(i+l)-Y'’’^(i-hl)=YM(i+l)-^(i)V(i)-0(i)U(i) 

Following  the  proposed  control  law,  the  control  effort  is  defined  to  assure  zero  value  of 
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the  predicted  error,  i.e. 

U(i) = 0(i)'[YM(i + l)-^(i)  V  (i)]  (2-34) 

Now  let  us  assess  the  "posterior"  error,  i.e.  the  system  error  which  takes  place  when 
control  effort  (2.34)  has  been  applied  at  the  i-th  step  of  the  control  algorithm. 

E(i+l)=YM(i+l)-Y(i+l)=YM(i+l)-RV(i)-QU(i),  or 
E(i+l)=YM(i+l)-RV(i)-Q©(i)  ‘[YM(i+l)-‘5>(i)V(i)],  or 
E(i+l)=YM(i+l)-Q0(i)-^YM(i+l)-RV(i)+Q©(i)-^^(i)V(i),  or 
E(i+l)=[I-Q©(i)']YM(i+l)-[R-Q®(i)'<^»(i)]V(i)  (234) 

Assume  that 

$(i)=R+A(i)  and  ©(i)=Q+S(i) 

where 

A(i)  and  S(i)  are  random  matrices  of  estimation  errors  and,  due  to  successful 
performance  of  the  RLSM  procedure,  their  norms,  l|A(i)||  and  ||5(i)ll,  are  effectively 
small.  Then 

Q©(i)-‘=Q[Q+S(i)r  and  I-Q©(i)'^  «0 

Similarly, 

R-Q©(i)  ‘$(i)  «  R-^(i)  =  R-R-A(i)  =-A(i) 

This  analysis  indicates  that  accurate  estimation/tracking  of  parameters  of  the 
controlled  plant  guarantees  the  feasibility  of  accurate  model-following. 
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3.  IMPLEMENTATION  AND  TESTING  OF  MR  CONTROL 
OF  OPTICAL  MIRRORS 

Implementation  of  the  developed  MR  control  systems  in  the  laboratory  conditions 
is  crucial  for  the  evaluation  of  the  feasibility  of  the  developed  technology.  It  allows  to 
obtain  convincing  information  on  the  advantages  and  disadvantages  of  particular  MR 
schemes  in  comparison  with  the  existing  control  methods.  The  following  claims  based 
upon  simulation  studies  are  to  be  proved  (or  disproved)  by  laboratory  tests  described  in 
this  section, 

MR  approach  allows  for  extending  bandwidth  of  particular  control  channels  of 
beam  steeres 

-  MR  approach  allows  for  the  elimination  of  "bending  modes"  from  transient 
response  of  control  channels  of  a  beam  steerer 

-  MR  reference  approach  results  in  the  elimination  of  cross-coupling  between 
control  channels  of  a  beam  steerer 

-  MR  approach  results  in  the  reduction  of  jitter  effects  on  the  beam  positioning 
accuracy 

-  MR  approach  results  in  the  increased  robustness  of  a  beam  positioning  system 

-  MR  approach  can  be  implemented  at  relatively  low  cost. 

Developed  control  techniques  were  implemented  using  the  laser  communication 
pointing-acquisition-tracking  (PAT)  testbed  of  Rome  Laboratory  equipped  with  500  Hz 
and  5KHz  optical  miaors  and  PC-based  control  systems  and  interfaces.  A  PC  equipped 
with  data  monitoring  software  package  was  utilized  for  generating  test  signals  and 
acquisition  of  system  responses  along  with  standard  data  probes.  Various  MR  schemes 
were  implemented  using  analog  and  digital  circuitry,  and  subjected  to  three  basic 
experiments.  The  first  experiment  was  designed  to  test  the  system's  ability  to  modify 
dynamic  response  of  a  single  control  channel.  The  second  experiment  tested  the 
decoupling  effects.  The  third  experiment  was  aimed  at  the  evaluation  of  the  jitter 
reduction  effect. 
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3.1  High  Bandwidth  Steering  Mirror 

The  PAT  testbed  [2]  provides  the  test  environment  for  the  high  bandwidth  mirror. 
The  hardware  of  the  PAT  system  includes  a  Hewlett  Packard  model-330  computer,  an 
electronics  rack  housing  the  support  electronics  for  the  electro-optic  and  electro¬ 
mechanical  components,  and  a  pair  of  optical  isolation  tables  upon  which  the  components 
of  the  primary  and  remote  stations  are  mounted. 

The  primary  station  optic  and  electro-optic  systems  can  be  categorized  into  six 
groups:  lasers,  detectors,  beam  conditioners,  ray  optics,  filters  and  beam  directors.  The 
primary  station  is  equipped  with  two  semiconductor  lasers  having  optical  output  powers 
of  10  mW  and  100  mW,  and  center  wavelengths  of  830  nm  and  800  nm  respectively. 

There  are  three  types  of  optical  detectors  in  the  primary  station,  two  types  for  tracking 
purposes  and  one  for  communications.  The  tracking  detector  is  a  quadrant  silicon 
avalanche  photodiode,  the  acquisition  detectors  are  improved  tetra-lateral  position 
sensitive  devices  (PSD)  and  the  communication  detector  is  a  silicon  avalanche  photodiode. 
The  beam  conditioners  and  ray  optics  assure  the  required  characteristics  of  the 
communication  and  beacon  laser  beams  for  detection  and  transmission.  The  optical  filters 
are  used  to  separate  these  beams  so  that  each  can  be  processed  individually.  The  beam 
directors,  coarse  steering  mirrors  (CSMl  and  CSM2)  and  the  fast  steering  mirror  (FSM), 
are  used  to  direct  the  transmitted  beam  and  to  position  the  received  beam. 

The  primary  station  electric  and  electro-mechanic  systems  include:  the  computer, 
the  computer  interface,  the  fine  steering  subsystem,  the  coarse  steering  subsystem,  and  the 
wide  field  of  view  (WFOV)  subsystem.  The  computer  interface  consists  of  a  HP  general 
purpose  input/output  (GPIO)  board  and  interface  electronics.  Both  the  fast  steering 
subsystem  and  the  coarse  steering  subsystem  are  analog  closed-loop  control  circuits.  The 
WFOV  camera  is  synchronized,  monitored  and  controlled  by  the  computer.”  [1] 

The  configuration  of  the  remote  station  is  similar  to  that  of  the  primary  station 
vrithout  the  acquisition  and  tracking  functions.  Two  remote  station  reside  on  a  single  4x8’ 
optical  table.  Remote  station  1  is  fixed  while  remote  2  is  mounted  on  a  rml  and  can  be 
driven  over  a  4  foot  distance.  The  beacon  and  communication  lasers  have  output  powers 
of  10  mW  and  50  mW,  and  center  wavelengths  of  780  nm  and  810  nm  respectively. 
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Optically,  the  tracking  laser  from  remote  terminal  2  propagates  across  the 
laboratoiy  (15  m)  impinging  on  the  #2  coarse  steering  mirror.  The  beam  is  relayed 
through  a  40x  beam  expander  to  the  fine  steering  mirror.  This  expansion  takes  place  to 
accurately  simulate  a  link  distance  between  a  low  Earth  satellite  and  a  geostationary  one 
(roughly  30,000  km).  From  the  fine  steering  mirror  the  beam  is  relayed  via  dichroic 
mirrors  through  the  communication  beam  pick  off  and  a  10  nm  optical  filter  and  400  mm 
lens  to  the  quadrant  detector.  Only  the  tracking  laser  was  used  for  these  measurements. 

The  control  system  for  the  coarse  steering  mirror  (CSM)  system  is  operated  in  a 
manual  mode  adjusting  the  beam  so  that  it  is  placed  in  the  field  of  view  of  the  quadrant 
detector.  The  CSM  is  further  adjusted  to  provide  a  beam  intensity  at  the  quad  cell 
sufficient  to  exceed  the  threshold  required  for  the  operation  of  the  fine  track  control  loop. 
Once  the  CSM  position  is  set,  the  quad  cell  is  adjusted  such  that  the  beam  is  centered  on 
it.  Figure  3-1  and  Figure  3-2  below  feature  voltage  vs.  position  characteristics  of  the 
quadrant  detector  employed  as  the  beam  position  sensor.  Note  that  zero  angle  occurs 
with  the  spot  at  the  edge  of  the  detector  for  the  axis  being  measured  and  centered  for  the 
other  axis.  Figure  3-2  shows  an  expanded  view  of  the  linear  region  (measurement  span). 
This  shows  that  a  measurement  accuracy  of  0.25  p  radians  is  possible.  It  is  desired  to 
operate  the  system  such  that  the  beam  remains  below  ±400  mV  so  that  it  remains  in  a 
linear  operating  region. 


Figure  3-1  Voltage  versus  angle  for  PAT  quad  cell. 
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Implementation  of  the  Reference  Model  and  Compensator.  The  reference  model 
and  compensator  were  implemented  in  analog  circuitry  according  to  the  transfer  functions 
derived  in  the  previous  chapter.  For  the  inverting  operational  amplifier  configuration, 
shown  in  Figure  3-3,  the  transfer  function  is  given  by  (using  the  infinite  gain 
approximation): 


Cl 


Figure  3-3.  Reference  model. 

The  component  values  were  chosen  from  standard  parts  available  in  the  lab  producing  a 
model  with  the  following  transfer  function: 

3030 

m,,(s)  =  m22(s)=  ~~~ 

iiv  /  22V/  s+3030 

Similarly,  the  transfer  function  for  the  lead  elements  of  the  compensator,  shown  in 


Figure  3-4,  is  given  by: 
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V„  _  R,(R,C3S  +  1) 
V.  R4(R5C4S+1) 


C3  C4 


Figure  3-4.  Lead  network 

Again,  component  values  were  chosen  from  standard  parts  producing  a  lead  network  with 
the  following  transfer  function: 

6.2E-5S  +  1 
6.2E-6S  +  1 


These  values  will  yield  a  lead  network  with  a  gmn  value  of  10  dB  at  a  frequency  of  8.1 17 
kHz. 

Cl  C2  Rll 


Figure  3-5.  Circuit  Diagram  of  single  channel  of  reference  model  and  compensator. 
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In  the  circuit  shown  in  Figure  3-5,  Ul,  Rl,  R2,  and  Cl  make  up  the  reference 
model.  The  output  of  the  reference  model  is  inverted  with  respect  to  its  input  and  is 
added  to  the  position  signal  from  the  mirror  via  C2  and  R3.  The  resulting  error  signal  is 
applied  to  U2,  U3,  and  U4,  which  represent  the  three  lead  elements  of  the  compensator. 
The  output  of  the  compensator  is  then  added  to  the  reference  signal  and  applied  to  the 
mirror  drive.  The  compensator  and  model  for  both  elevation  and  azimuth  were  built  on 
proto-boards. 

A  square  wave  input  from  a  HP  8116  function  generator  with  an  amplitude  of 
±200  mV  was  applied  first  to  the  input  of  the  azimuth  channel  and  then  to  the  elevation 
channel  and  the  outputs  were  measured  using  a  Phillips  digital  oscilloscope.  Figure  3-6 
and  Figure  3-7  show  the  response  of  the  system  with  the  conventional  control.  It  should 
be  noted  that  while  the  conventional  controller  deals  exclusively  with  the  “main  dynamics” 
of  the  controlled  plant,  the  model  reference  controller  improves  the  “main  dynamics”  and, 
in  this  case,  reduces  the  cross-coupling  and  the  effects  of  jitter. 


Figure  3-6.  Step  response  of  conventional  control  -with  azimuth  channel  excited. 
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Figure  3-7.  Step  response  of  conventional  control  with  elevation  channel  excited. 


Figure  3-8  and  Figure  3-9  represent  the  step  response  of  the  azimuth  and  elevation 
channels  using  a  reference  model.  Note  that  while  perfect  model  following  is  not 
exhibited,  clearly  the  dynamic  characteristics  are  the  same  or  better  than  the  original 
system.  The  overshoot  has  been  significantly  reduced  along  with  a  reduction  in  the  cross¬ 
channel  coupling. 


Figure  3-8.  Step  response  of  model  reference  control,  azimuth  channel  excited. 
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Figure  3-9.  Step  response  of  model  reference  control,  elevation  channel  excited. 


The  system’s  response  to  jitter  was  measured  through  the  introduction  of  acoustic 
vibration.  The  vibration  was  produced  by  applying  a  400  mV  sinusoid  with  sweeping 
frequencies  from  10  to  1000  Hz  to  a  audio  amplifier  and  speaker.  The  acoustic  vibration 
is  additive  to  the  jitter  due  to  atmospheric  transmission.  Figure  3-10  depicts  the  jitter 
spectrum  of  the  elevation  axis  due  to  the  atmosphere  alone.  The  atmospheric  jitter  in  this 
axis  was  much  worse  than  that  of  the  azimuth  axis  and  is  the  result  of  the  air  flow  patterns 
in  the  laboratory.  For  terrestrial  applications  of  free  space  laser  communication,  this  jitter 
represents  the  fine  steering  problem  and  seems  to  have  the  same  power  spectral  density  as 
that  of  satellite  jitter  [2].  The  spectra  shown  are  the  average  of  20  traces.  It  can  be  seen 
that  the  conventional  controller  offers  a  reduction  in  jitter  at  frequencies  below  300  Hz. 
The  model  reference  control  adds  little  additional  improvement  below  300  HZ  but 
effectively  adds  an  additional  7  dB  improvement  at  higher  frequencies.  The  addition  of 
the  acoustic  jitter  excites  various  mechanical  modes  in  the  table  as  demonstrated  by  the 
peaks  in  Figure  0-14.  Again,  when  fine  tracking  is  implemented  a  large  reduction  in  jitter 
is  seen  at  frequencies  below  300  Hz  along  with  an  equalization  of  the  modes  above  400 
Hz.  Once  model  reference  control  is  added,  a  5  dB  improvement  can  be  seen.  The 
presence  of  modes  at  60,  120,  and  180  Hz  is  due  to  electrical  noise  pick-up  in  the 
protoboard  used  to  implement  the  compensator  circuit. 
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Figure  3-10.  Atmospheric  jitter  spectrum  no  control  effort  applied 


Figure  3-11.  Atmospheric  jitter  spectrum  using  standard  feedback  control. 


Figure  3-13.  Composite  atmospheric  jitter  spectrum. 


Figure  3-14.  Acoustic  jitter  induced  spectrum  no  control  effort  applied. 


Figure  3-15.  Acoustic  jitter  induced  spectrum  with  standard  feedback  control. 
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Figure  3-16.  Acoustic  jitter  induced  spectrum  with  model  reference  control. 


Figure  3-17.  Acoustic  jitter  induced  spectrum  with  model  reference  control. 


3.2  Low  Bandwidth  Mirror: 

The  model  reference  control  technique  was  also  applied  to  a  Physik  Instrument 
model  S-320.2  three  axis  tilting  platform  that  provides  an  overall  steering  angle  of  ±1  mr. 

The  experimental  setup  used  a  collimated  HeNe  laser  passed  through  a  variable 
aperture  to  provide  a  collimated  beam  with  a  diameter  smaller  than  that  of  the  mirror. 

This  beam  was  then  directed  to  the  steering  mirror  and  then  through  a  1500  mm  lens  to  a 
RCA  30927E  quadrant  detector .  Figure  3-  below  shows  the  quad  cell  output  voltage  Vaz 


and  Vel  as  a  function  of  angle.  Unlike  the  PAT  system,  the  laser  was  not  modulated  nor 
were  electrical  or  optical  bandpass  filters  used  on  the  output  of  the  quad  cell.  In  order  to 
reduce  the  effect  of  ambient  light  induced  noise  in  the  system,  the  laser  power  was  kept 
high  and  the  room  lights  were  turned  down. 


The  control  for  this  system  was  implemented  digitally  using  a  Macintosh  Quadra 
650  computer.  The  control  system  is  shown  in  Figure  3-19.  The  shadowed  boxes 
represent  functions  which  take  place  in  the  DSP  card.  As  before,  M(s)  is  the  reference 
model  and  is  given  by: 


M(s)  = 


fm  „  (s) 


0 


o1 


m„(s)  0_ 


One  can  see  that  the  reference  model  does  not  exhibit  cross-coupling  and  cannot  be 
affected  by  environmental  jitter.  The  particular  choice  of  transfer  function 


mii(s)=m22(s)= 


4000 


s  +  4000 


This  model  is  a  little  faster  than  the  model  used  for  the  high  bandwidth  mirror  since  one  of 
the  goals  of  this  research  is  to  demonstrate  improved  performance  from  commercially 
available  devices.  The  controller  provided  by  Physik  Instrument  has  a  bandwidth  of  500 
Hz  and  does  not  provide  for  any  decoupling  between  axis.  As  previously  discussed,  the 
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mirrors  step  response  was  poorly  damped  and  highly  coupled.  It  was  necessary  therefore 
to  introduce  several  sections  of  lag  in  the  forward  path  and  an  explicit  decoupler.  They 
are  represented  by  H(s).  C(s)  is  an  additional  lowpass  filter  introduced  to  insure  that 
frequencies  above  the  Nyquist  rate  were  eliminated  before  conversion. 


Command 


Figure  3-19.  Piezo  mirror  control  system 


Once  again  Q(s)  is  a  3  stage  lead  network  providing  gain  of  64  dB.  Both  Q(s)  and  M(s) 
were  implemented  as  single  stage  HR  filters.  The  reference  signal  was  synthesized  in  the 
dsp  and  represents  a  ±200  mv  square  wave  with  a  period  of  21  Hz.  Position  signals  were 
measured  at  the  output  of  the  quad  cell  prior  to  filtering  by  the  Kronhite  analog  filter 
(C(s)). 

Figure  3-20  below  shows  the  open  loop  response  of  the  system  when  only  the  lag 
filter  and  decoupler  are  in  place  in  response  to  excitation  of  the  azimuth  channel.  The 
system  exhibits  a  first  order  response  wth  a  settling  time  of  about  5  ms.  There  is  very 
little  evidence  of  coupling. 
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Figure  3-20.  Open  loop  response  to  a  step  applied  to  the  azimuth  axis. 


Figure  3-21 .  Aamuth  step  response  when  feedback  alone  is  applied. 


For  the  case  where  only  feedback  is  applied  (Figure  3-21),  we  note  an  improvement  in 
performance  but  fail  to  achieve  good  model  following. 

With  model  reference  control  (Figure  3-22)  good  model  following  is  achieved. 
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Figure  3-22.  Azimuth  step  response  with  model  reference  control  applied. 


Acoustic  jitter  was  introduced  to  the  table  by  sweeping  frequencies  from  10  to 
1000  Hz  via  an  audio  amplifier  and  speaker  system.  Measurements  of  the  beam  position 
spectra  were  made  using  a  HP3561 A  dynamic  signal  analyzer  for  both  the  azimuth  and 
elevation  channels.  Since  the  data  were  similar  we  will  only  present  the  results  of  the 
azimuth  channel.  Figure  3-23  demonstrates  the  response  of  the  system  when  only  the 
commercial  controller  is  in  place.  There  are  prominent  modes  in  the  table  at  90  and  340 

Hz. 


Figure  3-23.  Acoustic  induced  jitter  control  is  provided  by  the  commercial  controller. 
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Figure  3-24  illustrates  the  improvement  in  jitter  rejection  when  model  reference  control  is  applied. 
It  can  be  seen  that  there  is  a  significant  improvement  over  the  conventional  controller. 


Figure  3-24.  Acoustic  induced  jitter  model  reference  control. 
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Figure  3-25  Comparison  of  commercial  versus  model  reference  control. 


3.3  Summary 

A  simplified  model  folloAving  control  was  demonstrated  on  both  a  high  bandwidth  fine 
steering  mirror  and  a  commercially  available  piezo-driven  mirror.  Both  analog  and  digital 
implementations  were  demonstrated  .  The  performance  of  both  systems  showed  an  improved  step 
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response  consistent  with  previously  obtained  simulation  results.  The  most  dramatic  improvements 
were  seen  in  the  piezo  mirror.  Its  performance  was  enhanced  to  the  level  of  that  of  the  very 
expensive  high  bandwidth  mirror.  This  conclusively  demonstrates  that  improvements  in  the 
control  system  will  allow  the  use  of  piezo  mirrors  for  beam  stabilization  in  lasercom  systems. 

The  ability  to  "enforce"  dynamics  of  a  time-invariant  reference  model  on  a  beam  steering 
mirror  indicates  that  the  developed  approach  results  in  the  enhanced  robustness  of  the  resultant 
control  system. 

While  reference  models  are  completely  protected  from  any  environmental  effects,  such  as 
jitter,  the  jitter  performance  of  the  system  should  be  improved  by  application  of  MR  controllers. 
This  improvement  was  successfully  demonstrated  with  a  7  dB  decrease  in  the  error  spectrum. 

This  is  particularly  impressive  at  frequencies  above  200  Hz  since  previous  methods  required 
physical  modifications  to  the  mirror. 

It  was  shown  that  implementation  of  MR  control,  in  addition  to  conventional  feedback 
controllers,  results  in  the  increased  efficiency  of  beam  steerers  without  the  penalty  of  increased 
costs  and  payloads. 
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4.  POWER  CONTROL 

Modification  of  the  laser  beam  steering  control  results  in  the  improved  positioning 
accuracy.  This  allows  to  assure,  with  some  degree  of  confidence,  that  laser  beam  power  is 
concentrated  within  a  much  smaller  angle  thus  improving  the  signal-to-noise  ratio  at  the  receiver 
cite.  Therefore,  the  minimum  received  power  requirement,  corresponding  to  the  acceptable  rate 
of  communication  errors,  can  be  satisfied  with  a  smaller  overall  transmitted  laser  power.  Once 
the  laser  communication  system  is  designed,  and  placed  on  orbit,  variations  of  the  laser  beam 
positioning  accuracy  can  be  utilized  for  incremental  changes  to  the  transmitted  power.  Operation 
of  the  laser  source  at  lower  than  nominal  power  level  results  in  improved  reliability  of  electronics 
and  reduced  cycling  of  electric  batteries.  It  has  been  proposed  to  optimize  the  transmitted  laser 
power  on  the  basis  of  the  amount  of  the  residual  jitter  [13, 15]. 

4.1  General  Considerations 

Assume  that  the  laser  source  of  a  transmitting  satellite  and  the  light  sensitive  device  of  a 
receiving  satellite  are  ideally  aligned.  Any  change  in  the  orientation  of  the  transmitter  with  respect 
to  receiver  can  be  represented  using  two  orthogonal  position  errors,  ex  and  ey,  defined  as  angular 
deviations  from  the  original  position.  It  could  be  seen  that 

-7i/2  <  ex  <  7i/2  and  -idl  <  ey  <  7i/2  (4.  l) 

and  the  laser  power  can  be  characterized  by  a  symmetric,  distribution  'P(ex,ey)  with  a  maximum 
at  ex=ey=0.  The  full  power  of  the  laser  source  can  then  be  defined  as 

n 

1 

Wo  =  J  J'P(e,„ey)de,.dey  (4.2) 

K 

Introduce  a  generalized  angular  position  error,  assuming  radial  symmetry, 

e=(ex^+ey2)l/2  (4.3) 

Represent  the  power  distribution  as 

'F(ex,ey)=Wo'F(e),  (4.4) 

where 
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n 

2 

J^(e)de=l 

n 

Assume  that  at  a  particular  distance  between  the  transmitting  and  the  receiving  satellites 
(which  is  irrelevant  at  this  point)  some  threshold  value  of  intensity,  5,  can  be  defined,  such  that 
any  intensity  below  6  is  not  sufficient  for  reliable  operation  of  the  communication  system.  Define 
the  maximum  position  error  corresponding  to  this  threshold,  by  solving  the  following 

equation; 

Wo'F(eMAx)=5,  or  ^(eMAx)=5/Wo  (4.5) 

Since  'P(e)  is  a  bell-shaped  distribution,  it  can  be  seen  that  a  decrease  of  the  SAVq  value  results  in 

the  increase  of  the  solution  eMAx-  This  is  not  surprising:  increased  power  of  the  laser  source 
always  results  in  the  increased  tolerance  of  any  laser  communication  system  to  position  errors. 

Now  consider  the  probability  density  of  the  position  error  which  can  be  described  by  a 
normal  distribution  with  mean  |i  and  vanance  a  : 

P(e)=0(e,o) 

The  probability  of  a  failure  in  the  communication  system  caused  by  position  errors  can  be 
defined  now  as 

00 

Pf  =  2  j^(e,o-)de 


One  can  realize  that  due  to  the  shape  of  the  normal  distribution,  a  decrease  of  the  standard 
deviation  o  results  in  the  decrease  of  the  failure  probability  Pp. 

Summarizing  the  above  discussion,  probability  Pp  can  be  expressed  as  some  function  of 
the  power  of  the  laser  source  and  the  standard  deviation  of  the  position  error. 

Pp=f(Wo,a) 


It  could  be  seen  that  derivatives  of  the  probability  of  failure 


^(W„,a) 


<0  and 


(4.7) 


Consider  the  following  equation,  presenting  the  "balance"  between  the  power  of  the  laser 


source  and  the  standard  deviation  of  the  communication  error; 
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°  da 


(4.8) 


where 


AWq  and  Ao  are  increments  of  the  power  of  the  laser  source  and  the  standard  deviation  of 
position  error. 

Equation  (4.8)  holds  for  small  increments  AWq  and  Ao  and  has  the  potential  of  laser 


power  control  on  the  basis  of  on-line  estimation  of  changes  Ao(t)  caused  by  satellite  jitter.  The 
derivatives  utilized  in  (4.8)  can  be  precalculated  and  used  as  coefficients  of  a  linear  control  law. 

The  following  equation,  is  not  limited  to  the  magnitude  of  increments  and,  therefore,  is 
more  convenient  for  engineering  calculations: 

f(Wo+AWo,o)+f(Wo,a+Ao)=0  ( 4.9) 

Although  an  analytical  solution  of  (4.9  )  is  possible,  it  constitutes  a  formidable  analytical 
task.  A  numerical  solution  of  (4.9)  can  be  obtained  by  specification  of  the  particular  standard 
deviation  of  the  position  error,  a*=a+Aa,  and  application  of  a  direct  search  procedure  for  the 
minimization  of  the  norm  |f(Wo+AWo,a)+f(Wo,o*)|  with  respect  to  AWq. 

The  block  diagram  of  Figure  4-1  below  illustrates  the  suggested  numerical  solution 
procedure. 
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Figure  4-1  Block  diagram  of  the  power  control  procedure 
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4.2  Description  of  a  Direct  Detection  PPM  System 

The  expression  for  the  number  of  received  signal  photons  for  a  M-ary  pulse  position 
modulation  (PPM)  scheme  can  be  written  as  [21] : 

Ks(Pt.Gt’®t)  =  (  410  ) 

where 

iIt,  tir  are  efficiencies  of  the  transmitter  and  receiver  optics  Gt,  Gr,  are  the  transmitter  and 
receiver  antenna  gains, 

Lt(Gt,0t)  is  the  pointing  loss  factor  associated  with  the  pointing  error  0t, 

T|  is  the  quantum  efficiency  of  the  detector, 
hv  is  the  photon  energy, 

Pt  is  the  average  power, 

T,  is  the  time  slot  period,  and 

it  is  assumed  that  the  peak  power  is  M  times  the  average  power. 

This  equation  can  be  recognized  as  very  similar  to  the  link  equation  introduced  in  chapter 
1  for  a  range  of  z  and  a  wavelength  X.  For  this  discussion,  consider  the  receiver  gain  based  on  an 
unobscured  circular  aperture  with  a  wide  enough  field  of  view  that  the  receiver  pointing  error  has 
a  negligible  effect  on  the  received  signal  power.  Given  these  assumptions  the  receiver  gain  is 
determined  by 


where  rR  is  the  radius  of  the  aperture. 

Similarly  for  a  transmitter  illuminated  by  a  Gaussian  beam  with  rms  width  W  and  an 
unobscured  aperture;  Gt  and  Lt  are  given  by: 


(4.12) 


Lt(Gt,0t)— exp(-GT0\).  (4.13  ) 

Figure  4-2  shows  the  number  of  photons  received  as  a  function  of  mispoint  angle  for  two 
power  levels.  The  minimum  signal  requirement  is  dictated  by  the  desired  probability  of  bit  error. 
The  intersection  of  the  minimum  signal  line  >vith  the  number  of  photons  vs.  angle  curve  defines 
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the  maximum  angle  before  a  burst  error  occurs.  It  is  not  surprising  that  the  higher  the  power  the 
larger  the  allowed  angle. 


Figure  4-2.  Number  of  photons  received  vs.  angle. 


The  two  axis  pointing  and  control  system  is  assumed  to  have  a  coarse  outer  loop  which 
will  respond  to  slow  changes  and  a  inner  fine  steering  loop  that  will  be  responsible  for  fast 
changes.  It  is  the  responsibility  of  the  coarse  loop  to  cancel  any  bias  in  the  system.  This  allows  us 
to  assume  that  the  single  axis  probability  for  mispoint  is  a  Gaussian  function  with  single-axis 
pointing  errors  0x  and  0y.  The  Gaussian  function  has  standard  deviation  Ox  and  Oy.  We  assume 

radial  symmetiy,  so  that  Oe  =  Ox  =  Oy  and  0  =  +0x  •  The  radial  distribution  in  normalized 

from  can  be  recognized  as  the  Rayleigh  distribution. 

P(0)  =  -^expf-;^  (414) 

Recall  that  by  definition  a  burst  error  occurs  when  the  signal  irradiance  falls  below  the 
level  required  to  maintain  the  average  probability  of  bit  error.  The  probability  of  burst  error  is 


given  by: 


f0  f-e'L. 
Pc  *  =  — r  exp  — ;  d0 


(4.15) 


where 
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0  =9*  represents  the  mispoint  angle  that  produces  the  minimum  number  of  photons 
necessary  to  close  the  link  for  a  given  probability  of  bit  error.  The  solution  of  equation  (4. 1 5  ) 
yields  the  reference  angle  in  terms  of  the  probability  of  burst  error  and  the  standard  deviation  of 
the  Gaussian  probability  distribution  for  the  single  -axis  pointing  error 

0*  =  aV2(-lnP^*)  (4.16) 

The  probability  of  error  for  an  M-ary  PPM  scheme  using  avalanche  photodiodes  as  the 
detector  is  given  by  [21,  22,  23]  : 


P.  = 


M 

1 

f®  1 

f  ^ 

o 

1 

M-l  T 

/4  V  > 

2(M-1) 

1“ 

1  2a;  J 

1  2of  J 

ay 

UX 

(4.17) 


where: 


2  2j.f  IbZk'l  j. 

Oo=GF^K3+—J-h  — 
of  =  G'f(^Ks+Kb+^J  + 


2kBT^Ts 


^  Ri.e' 


(4.18) 

(4.19) 


Po 


(4.20) 


p,  =g[ks+Kb+^)  +  ^  (4.21) 

The  parameters  G,  and  F  are  the  gain  and  the  excess  noise  factor  of  the  APD,  Ib,  and  I,  are  the 
gain  dependent  and  gdn  independent  APD  dark  current,  e  is  the  electron  charge,  ke  is  the 
Boltzmann’s  constant,  Teq,  Rl  are  the  detector  noise  equivalent  temperature  and  load  resistance 
respectively,  and  Kb  is  the  background  photon  count. 

The  solution  of  (4.17)  for  a  desired  average  probability  of  bit  error  and  given  background 
signal  yields  the  minimum  number  of  photons  necessary  to  close  the  link,  Ksmin-  This  value  can 
then  be  substituted  in  (4. 17)  which  can  then  be  solved  for  the  corresponding  maximum  mispoint 
angle  0*.  It  is  given  by  : 


0*  = 


G-p 


OtoGlPjy 


(4.22) 


where 
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is  a  constant  that  is  independent  of  the  transmitted 


power  and  transmitted  antenna  gain. 

Equating  (4. 16)  and  (4.22)  results  in  the  following  expression 

PT=^^expEGTa^(-lnPE*)]  (4.23) 

This  relationship  forms  the  basis  for  the  adjustment  of  the  transmitted  power  based  on  the 
measured  values  of  the  variance  of  the  pointing  error. 

In  this  case,  we  have  an  explicit  relationship  between  the  variance  of  the  jitter  and  the 
necessary  transmitted  power.  Therefore,  the  partial  derivatives  are  unnecessary  and  the 
transmitted  power  can  be  calculated  on  line  as  a  function  of  the  variance  of  the  jitter. 
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technologies  of  command,  control,  communications  and  intelligence  and  to 
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a.  Conducts  vigorous  research,  development  and  test  programs  in  all 
applicable  technologies: 

b.  Transitions  technology  to  current  and  future  systems  to  improve 
operational  capability,  readiness,  and  supportabiiity; 

c.  Provides  a  full  range  of  technical  support  to  Air  Force  Materiel 
Command  product  centers  and  other  Air  Force  organizations; 

d.  Promotes  transfer  of  technology  to  the  private  sector; 

e.  Maintains  leading  edge  technological  expertise  in  the  areas  of 
surveillance,  communications,  command  and  control,  intelligence,  reliability 
science,  electro-magnetic  technology,  photonics,  signal  processing,  and 
computational  science. 
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